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Q^ , We consider three active flavor neutrino oscihations where both the mass- 

. , square difi'erences play a role in atmospheric neutrino problem. We calculate 

00 . 

, the matter efl^ects arising due to propagation through earth. We demonstrate 

that these eff'ects improve the fit to the electron data vis-a-vis vacuum oscilla- 

, tions but make the fit to the muon data far worse, thus worsening the overall 

i fit. The results of our analytical calculation verify the numerical investigations 

I ' 

Oh! of this scheme presented earlier by Fogli et al. 
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I. INTRODUCTION 

There are three evidences for neutrino oscillations: 1) solar neutrino problem, 2) atmo- 
spheric neutrino problem and 3) the results of LSND experiment. Two flavor oscillation 
analysis of each individual problem gives three widely different values for neutrino mass- 
squared differences (A's). Oscillations between the three active flavor neutrinos have only 
two independent A's. Hence it was proposed that there should be at least one sterile neu- 
trino if one wants to explain all three evidences in terms of neutrino oscillations . However, 
various attempts were made to account for all three evidences in terms of three active flavor 
oscillations |@]. CHOOZ constraint on 13 mixing angle rules out most of these solutions 0. 

Recently three active flavor oscillations as explanation for all three evidences was revived 
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with the following scenario [^3- The larger A is about 0.3 eV^ which can explain LSND 
result. The smaller A is in the range 10~^ — 10^'^ eV^. The upper limit is set so that 
CHOOZ constraint on the relevant mixing angles will not apply to oscillations driven by 
this A. The lower limit is set by demanding that the solar neutrino survival probability 
should be independent of energy. In these three flavor oscillation schemes, the atmospheric 
neutrino oscillations are driven by both the A's. In particular, for the oscillations driven by 
the smaller A, the oscillation probability is almost zero for downward going neutrinos and 
is significant for upward going neutrinos. In particular these oscillations create the zenith 
angle dependence which is needed to explain the data. 

This scheme reproduces the overall suppression seen in the atmospheric neutrino data 
and also the zenith dependence of the ratio of muon to electron events. But it does not 
give a fiat dependence for the zenith distribution for electron events which is observed in 
the data. It also gives an incorrect zenith distribution for muon events. The large Ve ^ 
transitions driven by the smaller A, which are needed to explain solar neutrino problem, are 
responsible for this distortion of zenith angle distribution in atmospheric neutrino problem. 
In ref. [^] this scheme was investigated numerically. Their results are summarized in figs (4) 
and (5) of There it was also pointed out that matter effects worsen the already bad fit 
to data. Three salient features emerge from these graphs. 

• For both sub-GeV electron as well as muon events there is no appreciable difference 
in the event distributions with and without matter. This implies that matter and 
vacuum oscillation probabilities are the same. 

• For multi-GeV electron events matter greatly suppresses the large z/^ Ve oscillations 
and thus drastically reduces the large excesses of electron events which vacuum os- 
cillations generate especially at large distances of traversal, i.e for the upward going 
events. 

• Both for multi-GeV muon events and upgoing muon events inclusion of matter effects 
results in essentially a flat profile as a function of the zenith angle. 
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In this paper we show that these features can be explained by a simple analytic calculation 
of matter effects in this scheme by using perturbation theory. Very recently two scale 
oscillations of of atmospheric neutrinos were treated analytically [p|,p!0[ under a different set 
of neutrino masses and mixings. But these analyses can explain only LSND results and 
atmospheric neutrino data leaving solar neutrino problem untouched. 

II. TWO SCALE OSCILLATIONS IN MATTER- A PERTURBATIVE ANALYSIS 

We take Scheck-Barenboim (SB) scheme given in ||^ as a representative of models in 
which all evidences for neutrino oscillations are explained in terms of three active flavors. 
A recent update using latest reactor and accelerator data is given in [|n|]. The masses of 
the vacuum mass eigenstates are taken to be /ii, fi2 and fi^. The vacuum oscillation here is 
controlled by both the mass-squared differences A21 = yU^ ~ /^i A31 = fil — ftl- 
chosen to be about 0.3 eV^ to drive the oscillations seen at LSND and A21 is constrained to 
be a few times 10""^ eV^. Both scales average out to give an energy independent suppression 
in the solar neutrino case. The (12) mixing angle u should be close to 7r/4, so that the solar 
neutrino survival probability is 0.5. A recent analysis found that this energy independent 
solution of solar neutrino problem is allowed with good confidence level The (13) 

mixing angle is constrained to be quite small < 10° by reactor and LSND data. The (23) 
mixing angle ip is centered around 27 degrees to give a good fit to the ratio of muon to 
electron events in atmospheric neutrino data. 

The three flavor eigenstates are related to the three mass eigenstates in vacuum through 
a unitary transformation, 



where the superscript v on r.h.s. stands for vacuum. The 3x3 unitary matrix can be 
parametrized by three Euler angles {uj,(f),ilj) and a phase. The form of the unitary matrix 
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can therefore be written in general as, C/" = C^23('0) x Uphase x Uis{(f)) x Ui2{u;), where Uij{9ij) 
is the mixing matrix between ith and jth mass eigenstates with the mixing angle 9ij. The 
explicit form of U is 

/ \ 



JJV 



(2) 



where = sin0 and = cos0 etc. All the angles can take values between and 7r/2. [/" 
can also be written as 



Uel Ue2 Ue3 



JJV 



U^i U^2 Ufj^s 

^ Url Ut-2 Ut-3 j 

In the mass eigenbasis, the (mass)^ matrix is diagonal, 





A21 
A31 

In the flavour basis the (mass)^ matrix has the form 



V 



A31M31 + A21M21, 



(3) 



(4) 
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where 
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Matter effects can be included by adding A{r), to the e — e element of where 



A = 2V2(^) pE. 



(7) 



In the above equation m„ is the mass of the nucleon, the number of electrons per nucleon 
in the matter which is ~ |, and p is the density of matter in gm/cc. A can be written as 

A = 0.76 X 10"V X E. 

A is in eV"^, if E is expressed in GeV. In the atmospheric neutrino case most of the 
trajectories are through the mantle of the earth. We take the mantle density to be a 
constant of value 5 gm/cc. 

The matter corrected (mass)^ matrix in the flavour basis is 



where 



Mi = A31M31 + A21M21 + AM A, 



^ 1 0^ 



(8) 



Ma 



\ 



(9) 
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Now A31 ^ A21. Also for a typical neutrino energy in the multi-GeV range of about 3 
GeV A21 ~ A. Thus we work in an situation where A21, A ^ A31. 

In this approximation, to the zeroth order, both the matter term and the term pro- 
portional to A21 can be neglected in eq. (|). Then = A31M31, whose eigenvalues and 
eigenvectors are 



0;|i) 



0;|2) 



A31; |3) 
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We treat AMa + A21M21 as perturbation to the dominant term in and carry out de- 
generate perturbation theory. We first define two new states as follows 



|l') =a\l)+p\2) 
\2) = -(3\l) + a\2). 



(11) 
(12) 



a and /5 are determined by the conditions 
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a^ + /5^ = l, and (1 |// |2 ) = 0, 
where H' = AMa + A21M21. The condition eq. (0) leads to the following equation 



(13) 



a 



(14) 



1 - 2a2 (A21M) + {U!, - Ul,) 

Using the explicit form of the vacuum mixing matrix given in eqn. eq. ([I^) can be written 
in terms of mixing angles as 



a 



\ cos^ (j) sin 2uj 



1 - 2a2 ( A21 M) - (cos2 <j) cos 2uj) 
In the SB scheme uj is close to maximal mixing and is small so for this scheme we get 

A 



(15) 



a 



1 - 2a2 2A21 
From the above equation we get 



k{say). 



(16) 



a 



1 + 
1 - 



1 



(17) 
(18) 



For later use we also note that 
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The matter dependent mass eigenvalues are now given by 



rrin 



rrin 



A3i + Af/2, ~ A31. 



(19) 



Hence we get the mass squared differences in matter 

1 



Am 
^21 



A 



31 



A31, 



[A21 - A cos^ 



A 



32 



A 



31 



(20) 



By computing the first order correction to the old states one can derive the mixing angles in 
matter. It is straightforward to show that two of the mixing angles and ip are unaffected 
by matter (corrections of order {A/ and are negligible). This was also demonstrated in 
Pj. The (12) mixing angle uj is strongly modified by matter effects and is given by 



sm ujr. 



-(3 cos + a sin u; 



(21) 



Since u is constrained to be close to 45 degrees we get 



sm Ur, 



v/2 

1 
71 



{a -(3) 



1 



A 
A21 



(22) 



Let us estimate the quantity {A/ /S.2i)- Using the definition of A we get 

A 0.76 X 10"^ X 5 X E 



A 



21 



X X 10-4 



(23) 
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where we set p = 5 gm/cc, earth mantle density. The quantity x in the SB scheme typically 
is between 2 to 5. So one sees that (y4/A2i) — E where E is the neutrino energy in GeV. 
This implies 



Sm LJr, 



1 - 



E 



VTTe^. 

Now for typical multi-GeV events, E > 2GeV. This in turn implies 

E 



(24) 



1 



(25) 



From eq. (^) we see that sinum — > 0. In the case of two flavor mixing, this effect is known 
previously If the vacuum mixing angle is maximal, then the matter depedent mixing 
angle goes to zero at high energies. Here we have demonstrated that the same effect occurs 
in three flavor oscillations also, in the limit of perturbation theory being valid, sinum = 
implies = and hence the (12) scale in matter, A^, decouples from the oscillation 
probabilities involving electron neutrinos. Below we show that the same thing happens for 
electron anti-neutrinos also. The z/g oscillation probability in matter becomes 



PZ = '^U^^U^sUZU;:^,sin^ 1.27 



E 



sin^ 20 sin^ ip 



(26) 



as the rapidly oscillating scale A31 will average out to half. This transition probability 
is at most a few percent because of the smallness of 0. Let us now consider the case of 
antineutrinos. Since for antineutrinos A —A we get 



sm Ur 



1 + 



E 



VT+E^ 



In this case for typical multi-GeV energies we see from eq. (pT]) that smuj„ 
implies = 0. The z/^ z/g oscillation probability becomes 



(27) 
This 



— sin^ 20 sin^ ip 



E 



(28) 



which is again a few percent. 

Let us now contrast eq. (^6]) and eq. (|28| ) with the expression for the z/^ ^ z/g vacuum 



oscillation probability in the SB scheme. This is given by 

-2f/eif/e3t/^lt/^3 - 2U,2U,sU^2U^s, (29) 

where we have set the rapid oscillating term due to A31 equal to 1/2. In the SB scheme, 
c<j ~ 7r/4 and hence the magnitudes of Uei, Ue2, U^i and are 0.6 or more. Because of 
this, one gets large z/^ ^ transitions driven by A21, for large distances of travel (i.e. for 
upward going neutrinos). 

Let us now look at the survival probabilities for electron neutrinos and anti-neutrinos. 
By setting Um = for neutrinos and Um = 71/2 for anti-neutrinos, we obtain the matter 
dependent survival probabilities to be 

PJ^ = P^ = l-^sin220^L (30) 



We should constrast these with the vacuum survival probablilities 

P~ = P^e = 1 - 4f/eWe^2 (l-ST^) 

-2[/,2,f/f3 - 2U',,U',,. (31) 

The last two terms in the above equation are proportional sin^ (f) and are negligible. However 
the term containing A21 has a large coefficient and can cause significant oscillations for 
upward going events. 

Similarly, matter depdendent survival probability for muon neutrinos (anti-neutrinos) 
are obtained by setting Um = 0(7r/2). They turn out to be 

= 1 - 4 sin^ 2ip sin^ sin^ ( 1.27—^ J 

— 2 cos^ 0sin^ ■?/'(l — cos^ sin^ t/;) (32) 

The anti-neutrino survival probability has a similar expression, except that for anti- 
neutrinos is different from that of neutrinos. However, in both cases, the zenith angle 
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dependence is very weak because the oscillations driven by have a coefficient propor- 
tional to sin^ (p and hence are very small. In case of vacuum survival probabilities, there is 
siginificant zenith angle dependence coming from the oscillations of A21 whose coefficient 
^Uf,iU^2 is quite large. 

Using the oscillation probabilities derived above, let us calculate the zenith angle distri- 
butions of electron and muon events in atmospheric neutrinos. They are given by 

Ne = ^^P^e + ^ePee + ^-^Pj^e + ^e^ee, (33) 

where the products are to be understood as convolutions over energy. To include matter 
effects, the matter dependent probabilities should be used in the above equations. The zenith 
angle dependence of vacuum probabilities, coming from the oscillations driven by A21, gives 
the required zenith angle dependence to upward going muon events for multi-GeV energies. 
For downward going muons, vacuum oscillations in SB scheme predict larger suppression 
compared to what is observed. In case of electron events, the vaccum oscillation prediction 
for downward going events is in agreement with data but for upward going events, a huge 
excess is predicted, which again occurs due to large z/^ ^ oscillations driven by A21. When 
matter effects are included, we have seen that A^ either decouples, as in Pee, -Pgg, -P^ie and 
P^g or has a very small coefficient, as in P^^ and P^^. Hence, with the inclusion of the matter 
effects, the zenith angle distribution for both electron and muon event distrutions become 
very fiat. This is desirable from the point view of electron data but the flat distribution 
strongly disagrees with the muon data. Since the error bars in muon data are much smaller 
than those in electron data, getting a better ffi to electron data at the expense of muon 
data, worsens the overall ffi of the data. Thus we see that inclusion of matter effects in SB 
scheme worsen what was already not a good ffi. 

We see from fig. (3) of ref. that the matter corrected oscillations of SB scheme predict 
a mild zenith angle distribution for the upward going multi-GeV electron events. Whereas 
the upward going muon events are predicted to have a fiat distribution. These two features 
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can also be explained in terms of our calculation. For electron events, the visible energy is 
the same as the neutrino energy (> 1.33 GeV). For energy range of 1.33 — 2 GeV, at the 
beginning of multi-GeV events, uJm 7^ for neutrinos (and Um 7^ 7r/2 for anti- neutrinos). 
Hence in this energy range, A'^ does not quite decouple. Since the flux is higher at lower 
energy, these neutrinos lead to more events and hence we see some zenith angle dependence 
of upgoing multi-GeV electron events. In case of muons, for a visible energy of 1.33 GeV, the 
neutrino energy is greater than 1.5 GeV and hence the decoupling of is better because 
of the higher neutrino energy threshold. 

For the through going muon events the neutrino energy is much higher than typical 
multi-GeV energies. So again 

, ^ - 1 (34) 

and the discussions of the previous section apply and we will get a fiat distribution for the 
muon events when matter effects are included, which is seen in the last panel of fig. (4) of 

i 

Lastly we discuss the sub-GeV events. For these events A < A21. The parameter 
k = (y4/2A2i) is about 15%. For most of the energy range, the matter effects are negligible 
and one has essentially vacuum oscillations. This can also be seen in the first two panels of 
fig. (4) of ref. [g. 



III. SUMMARY AND DISCUSSION 

We did an analytical calculation of matter effects in a scheme where all the three evidences 
for neutrino oscillations are explained in terms of three active flavor oscillations. Our simple 
calculation reproduced most of the features of the matter effects in this scheme, which 
were numerically investigated earlier. Matter effects bring the predictions of this scheme for 
electron events closer to the data and the fit to the electron data is improved. But they wipe 
out the zenith angle dependence of multi-GeV muon events, which is in contradiction to the 
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data. Since the error bars of muon data are much smaller than those of electron data, the 
deviation away from the data makes the overall fit, with matter effects, to the data worse 
compared to the already bad fit of vacuum oscillation predictions of this scheme. 

Finally we wish to make a rather interesting comment. If in future these kind of schemes 
become favored phenomenologically, the main motivation for buildging a muon storage ring 
based neutrino factories is lost. At such factories, the neutrino energies are in the range 
10 — 50 GeV. We have seen that at such large energies, due to earth matter effects, Um — 
for neutrinos and Um — tt/2 for anit-neutrinos. In either case, the sub-dominant scale 
decouples from oscillations and CP Violation in lepton sector will be unmeasurably small 
in a neutrino factory environment. In such a case, the best option to observe CP violation 
then will be a conventional superbeam with a short baseline (of about 100 km) . 
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